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We perform the comprehensive analysis of the pentaquark (5Q) in / = 3/2 ± channel using 
anisotropic quenched lattice QCD. We employ the standard Wilson gauge action at j3 = 5.75 
and the 0(a) improved Wilson (clover) quark action on a 12 3 x 96 lattice with the renormalized 
anisotropy as a s /a t = 4. A large number of gauge configurations as iV con f = 1000 is analyzed, 
which is found to be essential to achieve a reliable measurement. We study three types of the 
Rarita-Schwinger intepolationg fields with 1 = 0: (a) the NK*-type, (b) the (color-)twisted NK*- 
type, (c) a diquark-type. As a result, we find only massive states as m$Q ~ 2.1 — 2.2 GeV in 
J p = 3/2~ channel, and m$Q = 2.4 — 2.6 GeV in J p = 3/2 + channel in the chiral limit. The 
analysis with the hybrid boundary condition (HBC) is performed to distinguish whether these 
resonances are compact 5Q resonances or two-particle scattering states. No low-lying compact 
5Q resonance states are found below 2.1GeV. 
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'The Monte Carlo simulations have been performed on the NEC SX-5 supercomputer at Osaka University. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos.sissa.it/ 



Anisotropic lattice QCD study of pentaquark baryons in spin 3/2 channel 



Takumi Doi 



1. Introduction 



The announcement of the discovery of the new-particle & + (1540) with narrow width by the 
LEPS group at SPring-8[[T]] has triggered enormous studies to understand this mysterious state. 
In fact, + has baryon number B = +1, charge Q = +1 and strangeness 5 = +1, and therefore 
has the minimal configuration of uudds, which is manifestly exotic. Although the impact of this 
discovery is undoubted, the experimental status is quite controversial: several groups confirmed the 
existence of while others reported the null results^. Note that even positive results have 

not determined the quantum numbers such as spin and parity. Numerous theoretical studies have 
been performed as well to understand the nature of + [[|]. In particular, one of the central issues is 
to realize the mechanism for the extremely narrow width as T <J lMeV. Among several scenarios 
proposed so far, J p = 3/2~ possibility^] is interesting to investigate. In fact, if ® + is J p = 3/2~ 
state, the decay to the KN scattering state becomes the d-wave. Because it is expected that the 
special configuration (0s) 5 is dominant in the ground-state in J p = 3/2~ channel, the decay to the d- 
wave is suppressed by the wave function factor. Note that there is further suppression by the d-wave 
centrifugal barrier, leading to the significantly narrow decay width. However, this mechanism has 
a possible disadvantage that such a state tends to be massive due to the color-magnetic interaction 
in the constituent quark models. J p = 3/2+ possibility attracts another interest considering the 
diquark model[]|]. In this model, J p = l/2 + is predicted for @ + , where + is composed of an 
anti s quark and two scalar diquarks with relative p-wave angular momentum. Narrow width is 
expected from this exotic structure as well as the p-wave centrifugal barrier. If QCD really exhibits 
this scenario, we expect that J p = 3/2+ state also exists as a LS-partner of & + . Under these 
circumstances, we present anisotropic lattice QCD results on 5Q states in J p = 3/2 ± channels. 
Note that while there have been several lattice QCD calculations of 5Q states[§, [7[ §, jfij [TT], p), 
13, 14, O. pj, these studies are restricted to J p = l/2 ± channels except for a very recent one fll6|]. 



2. Formalism 



We study the pentaquark (5Q) state analyzing the Euclidean two-point correlator as G^ v (t) = 
^j^i/a jU (t,x)i// v (0,0)^) , where i/fy is a Rarita-Schwinger interpolating field. In principle, we can 
use any field l/fy with appropriate quantum numbers. In practice, however, some fields may cou- 
ple strongly to a genuine 5Q state, while other fields may couple weakly. Unfortunately, there 
is no empirical knowledge about the structure of the 5Q state, i.e., about the suitable interpolat- 
ing field, we perform the comprehensive study using the three types of fields and examine how 
the results depend on the choice of the fields. We employ following isoscalar fields: (a) NIC- 
type, \fffi = £ abc {uJiCysdh) u c ■ {sdjfidd) + (m <-> d) (b) (color-)twisted NK*-type, which is an ex- 
tention to the one in Ref.[||], l/fy = e a b c {u^Cysdy) u c \ ■ (sdY^d c ) + (u <-> d) (c) diquark- type [[7]], 
which is an extention to the one in Ref . Jp7|] , yfy = £ a b c £def£cfg {u^Cjsdb) {u T d Cy^de) Cy=,s g . To 
analyze spin 3/2 state, we decompose the correlator as G;y(r) = P^ 2 ^G^ 3 ^(r) +p| / 1 ^G( 1 / 2 )(t), 
where G^I 2 \t) and G^l 2 \l) denote the spin 3/2 and 1/2 contributions to G(x), respectively, and 

(3/2) (1/2) 

P\j = 8jj — (l/3)YiYj,P)j = {\/3)jiYj. In order to study both parity states, the parity projec- 
tion is performed as well with P± = (1 =F Ya)/2. Note that all of the fields transform as i//,(t,x) — > 
— 74^,(t, — x)(i = 1,2,3) under the spatial reflection of quark fields as q(T,x) — > +74<?(t, — x). 
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In the 5Q state study from lattice QCD, one of the most important issues is the discrim- 
ination between two-particle scattering states and genuine 5Q compact resonances. In fact, in 
7 = 3/2, 7 = channel, NK and NK* scattering states are expected to have large contribution in the 
5Q correlator. To resolve this problem, we use two distinct spatial boundary conditions(BC), i.e., 
the periodic BC (PBC) and the hybrid BC (HBC), which is recently proposed in Ref.||T(}. In PBC, 
we impose the spatially periodic BC on u,d and s-quarks. As a result, all the hadrons are subject to 
the periodic BC, and their momenta are quantized as pi = 2n,-7r/L, (n,- € Z), where L denotes the 
spatial extent of the lattice. On the other hand, in HBC, we impose the spatially anti-periodic BC 
on u and d-quarks, whereas the spatially periodic BC is imposed on s-quark. Since N(uud,udd), 
K(us, ds) and K*(us, ds) contain odd numbers of u and d quarks, they are subject to the anti -periodic 
BC and the momenta are quantized as pi = (2n, + \)n/L. The crucial point is that the energy of 

the low-lying two-particle scattering state, E m in — yj m% + |j5 m in| 2 + J >n\ + \Pmin\ 2 (K = K or K*), 
is different between PBC and HBC because of the difference of the quantization of momenta. In 
fact, a drastic change is expected in the s-wave NK* state. In PBC, the lowest energy of NK* state 
is £ m f n c — niN + nix* . In contrast, since both N and K* are required to have non- vanishing mo- 
menta (±k/L, ±k/L, ±71 /L) in HBC, the lowest energy is raised up as £"™ c — ^mj f ±3n 2 /L 2 + 

\Jm\» ±3n 2 /L 2 , where the shift of energy amounts typically to a few hundred MeV for L ~ 2 
fm. The energy shift appears in NK(d-wave), NK(p-wave), NK*(p-wave) as well. Note that the 
minimum momentum in PBC is not p = 0, but p= (±2tf/L,0,0), (0, ±2n/L,0), (0,0, ±2%/L), 
because p = state can couple only to s-wave scattering states. Therefore, the lowest ener- 
gies are expressed as £™ n c ~ ^m 2 N ±An 2 /L 2 + ^m 2 R ±4K 2 /L 2 and £"™ c ~ ^Jm 2 N ±3x 2 /L 2 ± 

■J~rr?£ ± 3ti 2 /L 2 . These energy difference between PBC and HBC of scattering states contrast 
strongly with one in a compact 5Q resonance. In fact, because & + (uudds) contains even number of 
u and d quarks, it is subject to the spatially periodic BC not only in PBC but also in HBC, which al- 
lows & + to have p = in both of PBC and HBC. Therefore, the energy difference of a compact 5Q 
state is expected to be marginal, since it originates only from the change in its intrinsic structure. 
In this way, we can identify whether the observed states on the lattice are compact 5Q resonances 
or scattering states, by examining the energy difference between PBC and HBC. 



3. Lattice QCD parameters and Numerical results 

We calculate 5Q correlators using anisotropic quenched lattice QCD. To generate gauge field 
configurations, we adopt the standard Wilson gauge action at /3 = 5.75 on the 12 3 x 96 lattice 
with the renormalized anisotropy a s ja, = 4. Note that the anisotropic lattice is suitable for high- 
precision measurements of temporal correlators. The lattice unit is determined as aj l = 1.100(6) 
GeV from the Sommer parameter 1 = 395 MeV, and thus the physical lattice volume amounts to 
(2.15fm) 3 x (4.30fm). We use totally 1000 gauge field configurations to achieve the high statistics 
analysis. This is quite essential for our study, because the 5Q correlators for spin 3/2 states are 
found to be rather noisy. For the quark fields, we adopt the O (a) -improved Wilson (clover) action 
on the anisotropic lattice to suppress the discretization error. We adopt four values of the hopping 
parameters as K = 0.1210(0.0010)0.1240, which corresponds to m n jm 9 = 0.81,0.78,0.73 and 
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0.66 and roughly covers the region m s <, m <, 2m s . By keeping K s = 0.1240 fixed for s quark, we 
change K = 0.1210 — 0.1240 for u and d quarks for chiral extrapolation. In order to enhance the 
low-lying spectra, we use a gaussian smeared source (gaussian size of p ~ 0.4 fm) and point sink 



correlator in the Coulomb gauge. For details of the lattice QCD calculations, see Ref . [ |18| ] 



3.1 J p = 3/2 5Q spectrum 

We consider 5Q spectrum in J p = 3/2~ channel in PBC. In the upper side of Figs. |T](a), (b), 
(c), we show the effective mass plots for three fields, (a) the NK*-type, (b) the twisted NK*-type, 
(c) a diquark-type, respectively. Unless otherwise indicated, we use (k s , k) = (0.1240,0.1220) as a 
typical set of hopping parameters hereafter. In upper Figs. |l] (a) and (b), we find the plateaus around 
25 <, T <> 35 where the contamination from excited states as well as the backward propagation are 
suppressed. We perform a single-exponential fit in this region and obtain (a) /M5Q = 2.90(2) GeV, 
(b) wi5Q = 2.89(1) GeV, respectively. On the other hand, we see that the statistical error is too 
large in upper Fig. [T] (c), and therefore we do not use diquark-type field in this channel. A possible 
reason for such a large noise is that this field does not survive in the non-relativistic limit. Now, 
we perform the chiral extrapolation. Fig. ^| shows the 5Q masses in J p = 3/2~ channel against 
m n . Since the data behave almost linearly in mr % , we adopt the linear chiral extrapolation in rrr % . 
As a result, we obtain only massive 5Q states as (a) m^Q = 2.17(4) GeV, (b) m5Q = 2.11(4) GeV, 
respectively, which is too heavy to be identified with + (154O). Note that none of 5Q states appear 
below the NK threshold, although this threshold is raised up by about 200 — 250 MeV due to the 
finite extent of the spatial lattice, and the 5Q signal is expected to appear below the (raised) NK 
threshold considering the empirical mass difference between N+K(1440) and + (154O). 

To clarify whether our 5Q states are compact resonances or not, we perform the analysis with 
HBC. With the typical set of hopping parameters, it is expected that the s-wave NK* state is raised 
up by ~ 180 MeV, the d-wave NK state is lowered down by ~ 70 MeV, and a compact 5Q state 
shows no difference by switching PBC to HBC. Actual lattice simulations show m^Q = 2.98(1) 
GeV for both of NK*-type and twisted NK*-type fields, which corresponds to upper shift by 80- 
90MeV. Although the shift of m^Q is rather small, the value of is consistent with the s-wave 
NK* state within the statistical error. Therefore, we regard this state as an s-wave NK* state. In this 
way, we do not observe any compact 5Q resonances in J p = 3/2~ channel below the raised s-wave 

NK* threshold, i.e., in the region of E <, \Jm 2 N + p^ m + ^Jm 2 K , + p^ m , with \p m m\ — 499 MeV. 

3.2 J p = 3/2+ 5Q spectrum 

We perform the analysis for 5Q spectrum in J p = 3/2+ channel. The lower side of Figs. [I] (a), 
(b), (c), show the effective mass plots for (a) the NK*-type, (b) the twisted NK*-type, (c) a diquark- 
type, respectively. In all figures, we observe the plateaus and perform the single-exponential fit. 
We obtain (a) m 5Q = 3.34(3) GeV, (b) m 5Q = 3.11(4) GeV and (c) m 5Q = 3.16(2) GeV, respec- 
tively.Note that (a) agrees with the s-wave N*K* threshold while the latter two is very close to NK* 
threshold. Now, we perform the chiral extrapolation. In Fig. Q m^Q is plotted against m 2 n . Using 
the linear chiral extrapolation in m|, we obtain (a) m^Q = 2.64(7) GeV, (b) m^Q = 2.48(10) GeV 
and (c) m5Q = 2.42(6) GeV, respectively. Note that we observe again that all of m^Q data appear 
above the NK* p-wave threshold, which is located above the (raised) NK p-wave threshold. 
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(a) NK* (J =3/2") [PBC] 



<b) twisted NK* (J =3/2") [PBC] 



(c) diquarMJ =3/2 ) [PBC] 
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(b) twisted NK*(J p =3/2 + ) [PBC] 
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(c) diquark ( J P =3/2 + ) [PBC] 
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Figure 1: The 5Q effective mass plots in J p = 3/2~ (upper) and J p = 3/2 + (lower) channel in the periodic 
BC(PBC) for three types of fields: (a) the NK*-type, (b) the twisted NK*-type, and (c) the diquark-type. The 
dotted lines indicate the s(p)-wave NK* and the d(p)-wave NK threshold for J p = 3/2~(3/2 + ). The solid 
lines denote the results of the single-exponential fit performed in each plateau region. 

Next, we perform the HBC analysis. For J p = 3/2 + channel, the energy shift is somewhat 
minor change. With the typical set of hopping parameters, it is expected that the s-wave N*K* 
state is raised up by ~ 170 MeV, the p-wave NK* state is lowered down by ~ 60 MeV, the p- 
wave NK state is lowered down by ~ 70 MeV and a compact 5Q state shows no difference by 
switching PBC to HBC. As a lattice result, we obtain m 5Q = 3.38(2) GeV from NK*-type field. 
Although corresponding shift of 40 MeV is rather small, m$Q is again almost consistent with the 
s-wave N*K* threshold. Considering its rather large statistical error, this 5Q state is likely to be an 
s-wave N*K* state. The results from other two fields are very similar to each other. We obtain (b) 
m 5Q = 3.02(3) GeV and (c) m^Q = 3.08(4) GeV, respectively. These results correspond to lower 
shift by 80-90MeV, which is considered to be consistent with the NK* p-wave behavior. Therefore, 
these states are likely to be an NK* p-wave state. In this way, we observe no signal for compact 
5Q resonances. Although the assignments for plateaus are still afflicted by considerable size of 
statistical error, we can at least state that these 5Q states are all massive, which locate above NK* 
p-wave threshold and that these states are too heavy to be identified with ©+(1540). 



4. Summary and Conclusions 

We have studied J p = 3/2 ± pentaquark(5Q) baryons in anisotropic lattice QCD at the quenched 
level. We have employed the standard Wilson gauge action at j8 = 5.75 on 12 3 x 96 lattice with 
the renormalized anisotropy a s /a t = 4. The large statistics as N con f = 1000 has played a key role in 
achieving a solid result in our calculation. For the quark action, we have adopted O (a) -improved 
Wilson (clover) action with the hopping parameter as if = 0.1210(0.0010)0.1240, which roughly 
corresponds to m s <^ m <, 2m s . We have examined three types of the interpolating fields as (a) the 
NK*-type, (b) the twisted NK*-type, (c) the diquark-type. In J p = 3/2~ channel, we have ob- 
served plateaus in the effective mass plots except for the diquark-type field. Employing the linear 
chiral extrapolations in m|, we have obtained m^Q ~ 2.17 and 2.1 1 GeV in the chiral limit for the 
NK*-type and the twisted NK*-type correlators, respectively. We have performed the HBC analysis 
and have found that both of the observed 5Q states are s-wave NK* scattering states. In J p = 3/2 + 
channel, we have recognized plateaus in all the three effective mass plots. The chiral extrapolations 
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J P =3/2 J P =3/2 + 
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Figure 2: otsq in J p = 3/2~ (left) and 3/2 + (right) channel against nv\ for the three fields, i.e., (circle) the 
NK*-type, (box) the twisted NK*-type, and (triangle) a diquark-type. Data from diquark-type are not plotted 
in J p — 3/2~ because a plateau is not obtained in the effective mass plot. Open symbols denote the direct 
lattice results, while the closed symbols and the solid lines represent the results of the chiral extrapolations. 

have lead to m 5Q ~ 2.64,2.48,2.42 GeV for the NIC-type, twisted NIC-type and the diquark-type 
correlator, respectively. HBC analyses have been performed and the observed 5Q states in both 
of the twisted NIC-type and the diquark-type correlators are most likely to be NIC p-wave states. 
For the NIC-type field, it is most likely to be an s-wave N*K* state, although more statistics is 
needed to draw a definite conclusion. At any rate, whatever the real nature of these 5Q states may 
be, they are all considerably massive states in the physical quark mass region in both of J p = 3/2 ± 
channels, and cannot be identified as + (154O) without involving a significantly large chiral effect. 
For further studies, it is important to perform the systematic studies of the 5Q states using such as 

(I) unquenched full lattice QCD, (2) finer and larger volume lattice, (3) chiral fermion with small 
mass, (4) more sophisticated interpolating field reflecting the structure of 5Q, etc. 
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